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Abstract 

The decorrelation of the orientation of the event-plane angles in the initial state of relativistic Pb-Pb and p-Pb collisions, 
the “torque effect”, is studied in a model of entropy deposition in the longitudinal direction involving fluctuations of 
the longitudinal source profile on large scales. The radiation from a single wounded nucleon is asymmetric in space- 
time rapidity. It is assumed that the extent in rapidity of the region of deposited entropy is random. Fluctuations in 
the deposition of entropy from each source increase the event-plane decorrelation: for Pb-Pb collisions they improve 
the description of the data, while for p-Pb collisions the mechanism is absolutely essential to generate any sizable 
decorrelation. We also show that the experimental data for rank-four flow may be explained via folding of the elliptic 
flow. The results suggest the existence of long range fluctuations in the space-time distribution of entropy in the initial 
stages of relativistic nuclear collisions. 

Keywords: ultrarelativistic Pb-Pb and p-Pb collisions, event-by-event fluctuations, harmonic flow, event plane 
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1. Introduction 


During the collective expansion of the fireball formed 
in relativistic heavy-ion collisions azimuthal deformations 
of the density are transformed into azimuthal asymme¬ 
try of particle emission spectra mm- In the presence of 
collective flow, the particle spectra contain the harmonic 
components 


dN 

p±dp±dT] d(j) 


oc ...-kn2(pj_,77)cos[2((/i-V'2)] 

+ V3{P±,V) cos[3((/i - V'a)] + ■ • • 


( 1 ) 


In each collision, the event-plane of the second or third 
order harmonic flow is oriented predominantly along the 
direction of elliptic or triangular deformations of the fire¬ 
ball. It has been suggested that the angles of the event- 
plane orientation might vary as a function of pseudorapid¬ 
ity 0 or transverse momentum j3] . The effect leads to the 
factorization breaking for the two-particle cumulant flow 
coefficients, 

i4,A(qq2) < j (2) 


where U is the transverse momentum or pseudorapidity, 


VnA{tl,t 2 ) = ((C 0 s[n ((/>1 - </>2)])), (3) 
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Figure 1: Schematic view of the entropy distribution in an early 
stage of an ultrarelativistic nuclear collision. The matter deposited 
from each wounded nucleon occupies an interval in space-time ra¬ 
pidity with a randomly distributed end. As a result, the event-plane 
angles in the forward and backward bins are decorrelated. 


and the average is taken over events and over all particle 
pairs with particles i in a bin around ti. 

The factorization breaking in transverse momentum has 
been studied quantitatively in dynamical models mm in 
p-Pb and Pb-Pb collisions. The hydrodynamic response 
from fluctuating initial conditions can describe the experi¬ 
mentally observed event-plane fluctuations and the factor¬ 
ization breaking in pj_ mm- 

The decorrelation of the event-plane angles at differ¬ 
ent pseudorapidities is seen in a number of calculations, 
both in hydrodynamic, cascade, or hybrid models Uli- 
[H]. Nevertheless, a simultaneous description of the Pb-Pb 
and p-Pb data [S] poses a real challenge. In this paper we 
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propose a decorrelation mechanism which is capable to 
grab the basic experimental features of both reactions. A 
schematic view of the model is depicted in Fig. show¬ 
ing an early stage of the collision just after the two nuclei 
have passed through each other. The key ingredient is that 
the entropy deposition from the wounded nucleons m is 
made in string-like objects whose end-point is randomly 
distributed; some are longer and some shorter, with the 
length generated uniformly in the available rapidity inter¬ 
val. The idea is closely related to the model of Ref. [16] . 


2. The correlation measure 


It is very difficult to disentangle the genuine event- 
plane decorrelation due to the collective expansion of a 
“torqued” fireball from non-flow fluctuations of short range 
in pseudorapidity |3]. This difficulty is cleverly solved by 
using a factorization ratio using three bins with a large 
separation in pseudorapidity, as proposed by the CMS col¬ 
laboration |8]: 




VnA{-V\V^) 


(4) 


with the forward reference bin 4,4 < 77 ;, < 5 well separated 
from the two central bins where | 77 a| < 2.5. The departure 
of the factorization ratio from unity is a measure of 
the event-plane angle decorrelation as a function of the 
pseudorapidity separation Arj = 2r]a. 

In Ref. [11], the factorization ratio for elliptic and tri¬ 
angular flow in Pb-Pb collisions at ^/snn = 2.76 TeV is 
calculated in event-by-event viscous hydrodynamic sim¬ 
ulations with Glauber initial conditions. Assuming an 
asymmetric entropy deposition in space-time rapidity from 
left- and right-going wounded nucleons, one finds that the 
orientation of the fireball deformation depends on space- 
time rapidity, as the contribution to the fireball entropy 
from target and projectile wounded nucleons changes with 
rj [3]. Calculations show that the event-plane decorrelation 
in pseudo-rapidity can be described qualitatively, but the 
factorization ratio is noticeably underestimated. More¬ 
over, the calculation cannot reproduce the observed fac¬ 
torization breaking in p-Pb collisions. In the following, we 
discuss a mechanism introducing additional fluctuations in 
the entropy deposition, with long range correlations, that 
improves the description of the measured factorization ra¬ 
tio Tn- 

In the presence of collective expansion, the orientations 
of the event-planes and the elliptic or triangular deforma¬ 
tion are transformed into the orientation and the magni¬ 
tude of the corresponding harmonic flow components m- 
By the same mechanism, the torque of the event plane as a 
function of space-time rapidity is transformed into the ra¬ 
pidity dependence of the event-plane orientation extracted 
from particle spectra. This relation is expected to hold for 
decorrelation effects on large scales, while fluctuation in ra¬ 
pidity on small scales can be modified and washed out by 


the hydrodynamic evolution, resonance decays, mini-jets, 
etc. In the following, we investigate a model of fluctua¬ 
tions in the entropy deposition in space-time rapidity in 
the initial state. Hydrodynamic simulations show that the 
initial torque of the fireball in space-time rapidity is trans¬ 
formed into a very similar torque in the pseudorapidity 
dependence of the harmonic flow event-planes BM- Un¬ 
fortunately, precise hydrodynamic calculations including 
non-flow effects are very demanding. In this paper, event- 
plane decorrelation in pseudorapidity for the second and 
third harmonic flow are approximated by the event-plane 
decorrelation in spacetime rapidity in the initial state. 

Statistical hadronization, where a finite number of 
hadrons in a given bin is produced from a fireball with 
principal axes ipn, leads to large decorrelation effects [3] 
whose origin is trivial and needs to be canceled out. The 
CMS ratios Q accomplish this goal. Indeed, suppose we 
compute cumulants for the produced hadrons between the 
largely separated bins around rja and rji,. Then 

VnA{v‘^, v’’) = 


where the azimuths of the produced hadrons, and (f> 2 , 
are evaluated in some reference frame, ipniVa) and ipniVb) 
are the event-plane angles of the fireball, and 4>i and 
(j )2 are evaluated relative to tpniVa) and respec¬ 

tively. The factorization in Eq. ([^ applies if the torque 
angle magnitude is uncorrelated with the flow magnitude. 
The factors cancel out in appropriate ra¬ 

tios. For the symmetric A-A collisions the production 
around rja is the same as around —r]a, hence taking the 
ratio Q accomplishes the goal. For asymmetric colli¬ 
sions, as p-A, the appropriate measure proposed by CMS 
is yJrn{r]a,r]b)rn{-ria,-T]b). 

According to the above discussion, the factorization ra¬ 
tio can be written as 


, X ^ (cos[n(V>„(-??a) - ^njrib))]) 

Tn{Va,Tlb) (cos[n(7A„ (t?,) - V'n (%))]) 


( 6 ) 


where the average is taken over events. Expanding 
7/>„(±?7a) ^ '0n(O) ± yields 


rn{Va,1lb)- (7) 

(cos[n('0„(O) - V'n(??&))] -n sin[n(7/;„(0) - V’n(» 76 ))] Va) 
(cos[n(7/>„(0) -V'n (??&))] +n sin[n(7/^„(0) -'0„(?7(,))] 


For small values of the decorrelation angle, further expan¬ 
sion leads to 


rnir]a,r]b) - 1 - 2rt^((l/’n(0) -• (8) 

arj 

The deviation of the factorization ratio from 1 is found to 
be approximately linear in rja, as observed by the CMS 
collaboration [S]. The deviation of the factorization ratio 
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from 1 in the initial state is given by the correlation of the 
twist angle and its derivative 

l-r„ ~ 2n^((V'„(0) - ■i/’n(?7b)))-%^)?7a 

dr] 

Ot iii’niO) - 'ipn{Vb)f)Va ■ (9) 

The last proportionality holds approximately because of 
the strong correlation between '^^(O) — 4 ’niVb) and ■ 

The slope /„ of the linear dependence of 


rniVa, %) = 1 - ‘^fnVa (10) 

can be related to the variance of the event-plane angle 
difference between the central and the forward bin. Due to 
event-by-event fluctuations, ((V'n(O) — tpnivb)))"^) is found 
to be nonzero in several model calculations of the initial 
state niiiin]. 

The parameter used by the CMS collaboration, 

( 11 ) 


is approximately equal to the slope /„ of the linear depen¬ 
dence (101 for small factorization breaking. Parametrically 
oc rir, which amplifies the factorization breaking for 
higher harmonics. For centralities where the elliptic flow 
is strong one expects that the nonlinear contribution for 
dominates V 4 . In that case which leads to the 

relation 


asymmetric linear (averaged) distribution is used success¬ 
fully in the modeling of relativistic nuclear collisions [13- 
m]. Thus the model adopted by us to describe the fluctu¬ 
ations reproduces, upon averaging, the earlier approaches 
for observables computed from single rapidity bins. 

Fluctuations in the distribution of right- and left-going 
nucleons give a torque in the event-plane orientation even 
in the average model [3], that can partially reproduce 
the factorization ratios rnir]a,Vb) measured by the CMS 
collaboration. This average torque model predicts, as 
we shall see, a very small factorization breaking in p-Pb 
collisions, unlike observed experimentally. Let us note 
that similar effects are expected in string models p3] . 
if rapidities of the color charges at the ends of the flux 
tube fluctuate. The investigated mechanism is restricted 
to fluctuations which are long-range. The presence of 
any additional torque S-ipniv) of th® event-plane angles 
+ 5ipn{,'n)^ coming from local clusters, thermalized 
jet remnants, etc., would not modify the factorization ra¬ 
tio TniVa, Vb) (sse the discussion of the preceding Section), 
if the production in the forward and central bins is un¬ 
correlated, {StpniFVa)S'il^niVb)) = 0. The same argument 
applies for short-range non-flow correlations, as pointed 
out by the CMS collaboration |S]- 

The emergence of the torque effect relies on two features: 

1. asymmetric source profile in pseudorapidity, and 

2. fluctuations. 


Fi74 Fl (12) 

well satisfied in the experiment |S]. Relation 0 cannot 
be easily applied to compare the size of r 2 and as the 
correlation between 7n(0) — '^’ni'Hb) and is stronger 

for n = 2 than for n = 3. 

3. Torque model with strings of fluctuating length 

The entropy distribution in space-time rapidity in 
the initial stage is not yet fully understood. In 3-1-1- 
dimensional hydrodynamic calculations, the initial pro¬ 
file in the longitudinal direction is often assumed as a 
smooth symmetric function. This assumption is sufficient 
to obtain an average description of pseudorapidity spec¬ 
tra in symmetric collisions. However, the radiation from 
forward- and backward-going color charges naturally leads 
to asymmetric distributions in rapidity m- Following this 
idea, we assume a simple model where gluons radiated 
from a charge moving with rapidity yt, are distributed in 
rapidity uniformly in a range [ya, yb], with the end position 
ya taken as random (cf. Fig. [^. When the distribution of 
ya is uniform in the available range [-j/beam, J/beam], which 
is what we assume, then the averaged distribution has a 
linear dependence on rapidity. Notably, such an approxi¬ 
mately linear dependence of the density of particles emit¬ 
ted from a single wounded nucleon has been identified from 
particle spectra in asymmetric d-Au collisions [TS]. The 


The fluctuations included in our model are two-fold. First, 
we incorporate the discussed fluctuations of the emission 
profile, as depicted in Fig. second, we fluctuate the 
strength of the sources, overlaying a suitable distribution 
over the distribution of the wounded nucleons. The com¬ 
bined amount of fluctuations is controlled by the multi¬ 
plicity distributions. In particular, in p-Pb collisions we 
set the parameters of the overlaid distribution in such a 
way that we reproduce the CMS data in Fig. 

Our calculations are carried out with GLISSANDO |^[23] . 
The realistic NN inelastic collision profile for the LHC en¬ 
ergies is taken from Ref. [25]. We use an excluded distance 
d = 0.9 fm when generating the nucleon configurations in 
the nuclei. The total inelastic NN cross section is 64 mb for 
Pb-Pb collisions at ,/snn = 2.76 TeV and 70 mb for p-Pb 
collisions at y^SNN = 5.02 TeV. The fireball density in the 
transverse plane and pseudorapidity rj is taken in a form 
as a sum over the iV+ right-moving and N~ left-moving 
wounded nucleons, 

N+ N- 

s(x, y,r])=Y^ gf (x, y,g) + J2 9i(.X, y, rj). (13) 

1=1 i=l 

The source density in the transverse plane involves a super¬ 
position of strength; Wi is the superposed random weight, 
described in more detail below. We take into account the 
admixture of the binary collisions [351 [37] . If denotes 
the number of collisions of the z-th nucleon with^tjie nu¬ 
cleons from the other nucleus, then W™" = J 2 j=i is 
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Figure 2: Multiplicity distribution in p-Pb collisions, where the data 
are for charged tracks with > 0.4 GeV and \r]\ < 2.4 measured 
by CMS |30| . and the line denotes the corresponding results of the 
torque model with rapidity fluctuations, convoluted with a negative 
binomial distribution for the strength of the Glauber sources. 


the acquired random weight for the binary component. A 
necessary smearing is achieved with a smoothed Gaussian 
form centered around the position of the nucleon, (xi^yi). 
Combining these elements yields 


9tix,y,v) = [w,{l-a)h^{T]) + Wi°'^a]H{r]) 


(14) 


The width of the smearing Gaussian is tr = 0.4 fm, and 
the mixing parameter controlling the contribution of the 
binary collisions is a = 0.15. 

The longitudinal density profile, according to the earlier 
discussion, has the form 


hi,±{v) = 26»[±(?7* - 77)], 


(15) 


where 9 denotes the step function and iji is distributed 
randomly in the range [-j/beam, 2/beam]- Upon averaging 
over rji, the function (15) yields “triangular” distributions 


used successfully in previous works. The overall rapidity 
profile i?(?7||) is taken as a plateau with Gaussian tails |28) . 


H[rf) = exp ( — 


- ^?p)^Q(l^l - ^p) 


20-2 


(16) 


where 0 ^ = 1.4 and rjp = 2.5 PU] . 

To summarize the above construction, the wounded nu¬ 
cleons lead to fluctuating and asymmetric production in 
pseudorapidity, while the binary collisions yield symmet¬ 
ric emission. 

To set the parameters of the overlaid distribution pro¬ 
ducing the weights Wi we proceed as in Ref. [H] . The pro¬ 


duction of charged particles from each source of Eq. (14) 
is described by the negative binomial distribution 


^A.K(n) = 


r(r 


k)A”k« 


r(K)n!(A -I- «;)"+'' 


(17) 


where the hadron multiplicity n has the mean and variance 
given by A and A(1 -|-A/k), respectively. In Fig. [^we show 
the result of the model fit to the GMS data [30]) where 
a very reasonable agreement in the large multiplicity tail 


is obtained. The optimum parameters are A = 4.6 and 
K = 1.4. We note that without the fluctuations in rapid¬ 
ity K = 0.9 [30]) i-S., the variance of the fluctuations of the 
strength of the sources must be larger in this case to repro¬ 
duce the same distribution of hadrons. Assuming that the 
statistical hadronization following the deterministic hydro- 
dynamic phase brings in an additional Poisson distribution 
for the number of hadrons, the weights Wi of entropy of 
the Glauber sources follow the T distribution [28] . 


Pv{w) 


r(^) 


(18) 


While the description of the model presented in this sec¬ 
tion seems rather involved, we note that apart for the fluc¬ 
tuations of the longitudinal extent of the sources from the 
wounded nucleons, which is novel and which upon averag¬ 
ing yields the previously used emission profiles, the other 
elements (admixing binary collisions, overlaying the T or 
negative binomial distributions) are standard in state-of- 
the art modeling of the Glauber phase of the collision, and 
the model parameters are fixed in the same way as in pre¬ 
vious studies. 


4. Results 

The factorization ratios for the second and third har¬ 
monic in Pb-Pb collisions are shown in Fig. [^ The cal¬ 
culation in the torque model with long-range fluctuations 
in pseudorapidity describes very well the data for semi¬ 
central and peripheral collisions, both for and r-j,. We 
note, however, the lack of agreement for r-i in central col¬ 
lisions, where the decorrelation is significantly overesti¬ 
mated in the model. The result is nontrivial, as we were 
not able to adjust the emission profile used in the model 
to improve without spoiling the agreement for r^. The 
data of the GMS collaboration show that the factorization 
ratio in central collisions depends on the choice of the ref¬ 
erence pseudorapidity bin rji, [5]. This indicates that in 
central collisions fluctuations in the rapidity distribution 
of short range or non-flow correlations become relatively 
more important. Such correlations could originate from 
hard physics that is outside our model of the initial stage. 
This issue calls for further studies. 

The effect of fluctuations in the entropy distribution in 
rapidity is most striking for p-Pb collisions, as demon¬ 
strated in Fig. (|^. The experimental data show a signif¬ 
icant factorization breaking in the second harmonic. The 
calculation in the model with average entropy profile in 
rapidity (dashed-dotted line in Fig. gives almost no fac¬ 
torization breaking. The reason is simple, as the transverse 
profile at a given space-time rapidity is dominated by the 
contribution from wounded nucleons from the Pb nucleus. 
When all these sources deposit the entropy in space-time 
rapidity in a similar way, the event-plane orientation will 
show almost no rapidity dependence. The picture changes 
dramatically if the entropy density in space-time rapidity 
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Figure 3: Factorization ratios rn for the elliptic and triangular flow, 
plotted as functions of the central pseudorapidity bin position, ob¬ 
tained from the torque model with fluctuating entropy distribution in 
rapidity (lines), and the data of the CMS collaboration (symbols) [8]. 
In panels (a), (b) and (c) are presented results for centralities 0-5%, 
20-30% and 50-60%, respectively. 


for each source fluctuates, hence becomes different for each 
of them. Since the contribution of each source from the 
Pb nucleus to the entropy density at a given space-time 
rapidity varies, in consequence the event-plane orientation 
between the forward and central bin decorrelates notice¬ 
ably. The calculation in the torque model with fluctuating 
entropy distribution in rapidity describes surprisingly well 
the experimental data (dashed line in Fig. |^. This is the 
key result of our paper, which shows that the incorporation 
of the long-range pseudorapidity fluctuations is crucial to 
explain the large decorrelation seen in p-Pb collisions. 

In Fig. we show the parameter of Eq. 0 for dif¬ 
ferent centralities. The calculation describes properly the 
measured r 2 and in Pb-Pb collisions from semi-central 
to peripheral collisions. In central and ultra-central colli¬ 
sions the model overestimates the factorization breaking, 
especially for r 2 . As stated above, this indicates that in 
central collisions other sources of correlations appear that 
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Figure 4: Factorization ratio for the elliptic flow in p-Pb collisions, 
plotted as a function of the central pseudorapidity bin position. 
The data from the CMS collaboration (symbols) agree well with 
the torque model with fluctuating entropy density in pseudorapid¬ 
ity (dashed line), while the model with averaged entropy profile in 
rapidity (dashed-dotted line) fails spectacularly. 
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are not captured in our model. The calculation repro¬ 
duces r 2 measured in p-Pb collisions, but not its centrality 
dependence. The experimental data for the fourth order 
harmonic flow F^, scaled by a factor 1/4 are very close to 
This is in agreement with relation 
and is consistent with the collective flow scenario. 


the numbers for F?. 

(5, 


We stress that this relation is independent of the specific 
model of initial conditions, and holds only under the as¬ 
sumption that collective flow with large values of V 2 is 
generated. 


5. Conclusions 

In this paper the factorization breaking for event-plane 
angles defined at different pseudorapidities has been inves¬ 
tigated with the help of the factorization ratio coefficients 
r„ [S]. The observed factorization breaking confirms qual¬ 
itatively the existence of event-plane decorrelation, as sug¬ 
gested already in Ref. [3] . We have shown that the strength 
of the factorization breaking is a sensitive measure of the 
fluctuations of entropy deposition in space-time rapidity. 
The CMS data for p-Pb collisions strongly suggest the ex¬ 
istence of such fluctuations, that are independent for each 
wounded nucleon. We have studied a simple model of fluc¬ 
tuating entropy distribution in space-time rapidity, where 
the entropy production profile from wounded nucleons is 
approximately uniform, but the position of its end-point 
in pseudorapidity fluctuates. Calculations within the pro¬ 
posed torque model with fluctuating entropy distribution 
describe fairly well the measurements of r 2 and in Pb-Pb 
collisions, except for central collisions. The fluctuations in 
the pseudorapidity profile of the initial fireball are abso¬ 
lutely essential in reproducing the data for p-Pb collisions. 

We have shown on general grounds that the rank-4 slope 
coefficient is very close to 4 F 2 , which is confirmed by 
the data and which is one more signature of collectivity in 
the fireball evolution. 

We thus argue, based on our analysis, that the measure¬ 
ment of CMS collaboration of the factorization breaking 
for different pseudorapidities demonstrates the existence 
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Figure 5: Effective “slope” parameter of Eq. 0 plotted as a 
function of the number of charged tracks. The data from the CMS 
collaboration (symbols) [8] are compared to calculations in the torque 
model with pseudorapidity fluctuating entropy density (lines) for 
and F^. The data points for F^ are scaled by 1/4 (full triangles) to 
test the relation & 
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of fluctuations in the initial fireball density. A successful 
description of the observed collective flow requires the in¬ 
troduction of realistic fluctuating entropy distributions in 
the initial stage of the hydrodynamic evolution. The pro¬ 
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very interesting to have similar estimates from the color 
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